Lifetime reduction of surface states at Cu, Ag and Au(lll) caused by impurity 
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We present density-functional results on the lifetime of the (111) surface state of the noble metals. 
We consider scattering on the Fermi surface caused by impurity atoms belonging to the 3d and Asp 
series. The results are analyzed with respect to film thickness and with respect to separation of 
scattering into bulk or into surface states. While for impurities in the surface layer the overall 
trends are similar to the long-known bulk-state scattering, for adatom-induced scattering we find a 
surprising behavior with respect to the adatom atomic number. A plateau emerges in the scattering 
rate of the 3d adatoms, instead of a peak characteristic of the d resonance. Additionally, the 
scattering rate of Asp adatoms changes in a zig-zag pattern, contrary to a smooth parabolic increase 
following Linde's rule that is observed in bulk. We interpret these results in terms of the weaker 
charge-screening and of interference effects induced by the lowering of symmetry at the surface. 

PACS numbers: 73.20.At,73.20.Hb,73.50.Bk 
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I. INTRODUCTION 

The lifetime of the noble- metal (111) surface states 
has been the subject of numerous theoretical and ex- 
perimental investigations. These states are formed due 
to the breaking of translational symmetry at the sur- 
face in combination with the band gap existing in the 
surface-projected bandstructure in the [111] direction. 
Due to their localized nature, they play a role in a 
variety of surface-related phenomena, such as cataly- 
sis, adsorbate interactions, or surface transport, and can 
be probed by spectroscopic surface-sensitive techniques 
such as scanning tunneling microscopy and spectroscopy 
(STM and STS), photoemission, etc..— They furthermore 
constitute the fundament for the appearance of quantum- 
interference phenomena such as the famous quantum mi- 
rage in quantum-coral systems^ 

Although surface states are orthogonal to bulk states 
when only the periodic crystal and surface potential is 
considered, in reality there are always coupling mecha- 
nisms of the surface states to each other and to the bulk 
states. Mainly, these mechanisms are electron-electron 
interactions (including interactions with collective elec- 
tron modes such as plasmons), henceforth labeled "ee," 
electron-phonon interactions, henceforth labeled "ep," 
but also scattering at defects of the surface such as step 
edges, islands or adsorbate atoms, henceforth labeled 
"ed." The total lifetime r associated with these mecha- 
nisms is related to a sum of the scattering rates by 
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At low temperature, the phonon contribution becomes 
negligible, while the contribution of the electron-electron 
interaction depends on the electron energy E with respect 
to the Fermi level Ep .— For E close to Ep , and assuming 
Fermi-liquid behavior, the contribution to the decrease of 
the lifetime from electron-electron interactions behaves 



as t^ 1 ~ (E — Ep) 2 and is thus also ultimately frozen 
out. However, the scattering at defects remains finite. 

The lifetime r of surface states is an experimentally 
measurable quantity. Most commonly used methods for 
this purpose are angular resolved photoemission spec- 
troscopy (ARPES^ and STS ARPES can detect the 
full surface band structure, E^ and (with 1/rg appear- 
ing as a broadening T of the electron bands) , but defect 
scattering is necessarily averaged over a large surface area 
and cannot be disposed of. In fact, one of the first ex- 
periments where the importance of defect-scattering was 
detected was based on ARPES 4 It was shown that the 
lifetime r of the Cu(lll) surface state actually decreased 
for energies close to Ep, contrary to the predictions of 
Fermi-liquid theory; this effect was attributed to strong 
defect-scattering at those energies. 

STS, on the other hand, follows the onset of the sur- 
face state in the dl/dV signal which is proportional to 
the density of states, changing shape from step- like (as it 
would be for a two-dimensional electron system), to ob- 
taining a certain width T ~ r _1 due to scattering Ad- 
ditionally, r can be found via the phase coherence length 
L<ji by probing a quantum interference pattern in the 
dl/dV curves produced by nanoscale resonators*^— In 
particular, STS has the advantage that one can choose a 
defect-free surface region, in order to isolate the electron- 
electron and electron-phonon scattering from the defect- 
scatteringJ&M^ 

Whereas in most of the experiments localized defects 
have been avoided as much as possible, there are a few 
experimental studies which followed the idea of control- 
ling the defect concentration, in order to extrapolate to 
a defect-free surface. The first of these experiments has 
been performed with ARPES for different concentration 
of potassium atoms on Cu (111)^ where the structural 
disorder was investigated and quantified with LEED. In 
references fl7rll9l . similar experiments have been reported, 
in which the surface-state lifetimes at several coverages 
of Cu adatoms on Cu (111) have been measured in order 
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to extrapolate to zero-defect concentration. 

Theoretical approaches to the surface-state lifetime ad- 
dress, in most cases, the electron-electron and electron- 
phonon interaction, under the assumption of an ideally 
defect-free surface. Progress has been made in calcula- 
tion of the effects of the electron-electron interaction, as 
the initial degenerate electron gas model 2 ^ was replaced 
by the more elaborate GW approximation^! of many- 
body theory. Applied to image potential states^ this 
approach led to a nice agreement with experiment,— ac- 
counting for both intra-band and inter-band scattering. 
Concerning clectron-phonon scattering, already calcula- 
tions based on the simple Debye model lead to a sat- 
isfying agreement between numerical and experimental 
data^ while more accurate ab initio calculations of the 
phonon density of states comprise the investigation of 
bulk and surface phonons . 23 ' 24 Generally, at low temper- 
atures the electron-phonon contribution of the scattering 
rate behaves as T = I-kXIibT with A being the electron- 
phonon coupling parameter. This linear behavior is also 
found by experiment S£ 

However, in particular on scattering of surface states 
at adsorbate atoms, only few calculations have been pub- 
lished to date, mainly based on the wavepacket propaga- 
tion method, which is a real-space, real-time approach^ 
A systematic theoretical ab initio investigation of life- 
time reduction caused by (elastic) scattering at differ- 
ent species of adatoms or substitutional impurities in the 
(111) surfaces of the noble metals is still lacking. Bridg- 
ing this gap is the purpose of this paper. Furthermore, 
many of the aforementioned lifetime measurements (es- 
pecially in STS experiments) have been performed at 
the lower band edge of the surface-state band, where 
the contribution to the surface-state linewidth caused by 
defect scattering merges with the electron-electron and 
the electron-phonon contribution. This is different at 
the Fermi energy (at least at low temperatures), where 
defect scattering becomes the dominant mechanism. In 
addition, the Fermi level is the relevant energy for trans- 
port processes. These are the main reasons, why we have 
performed our calculations at the Fermi level. 

The paper is organized as follows: We start with 
an introduction to the theoretical formalism, i.e. mul- 
tiple scattering theory, based on the Korringa-Kohn- 
Rostoker (KKR) Green function method. We proceed 
with an analysis of computational and numerical aspects 
of the calculations, demonstrating differences between 
full potential calculations and those performed within the 
atomic sphere approximation. We study the lifetime re- 
duction due to scattering at 3d and Asp impurities in the 
first surface layer and in the fee adatom position on the 
noble metal thin films Cu, Ag and Au, which turn out to 
show similar trends. While the trend for the scattering 
rates for impurities in the first surface layer is found to be 
similar to that calculated in bulk materials 2 ^ (known also 
from residual resistivity results 2 ^), the trend for scatter- 
ing at adatoms remarkably differ. Therefore, we focus on 
adatoms, and analyze the trend of scattering rates for 3d 



and Asp adatoms on Cu(lll). This analysis includes an 
angular-momentum resolved study of surface-state life- 
times as well as a resolution of scattering into bulk and 
surface states. Finally, the trend observed for magnetic 
impurities and adatoms is investigated. 

A note should be made here on our approximation con- 
cerning magnetic impurities. Below the Kondo tempera- 
ture, the magnetic moment of these defects is quenched 
by many-body effects^ The Kondo temperature varies 
over orders of magnitude with respect to impurity type 
and host, depending exponentially on the position of the 
d resonance and on the hybridization. Therefore, there 
is no 'unique' temperature above which all 3d impurities 
become simultaneously magnetic. At a sufficiently low 
temperature, however, they should all be non-magnetic 
in the sense of the Kondo effect. In this case, the non- 
magnetic density of states in the local-density approxi- 
mation does not represent the physical density of states, 
except exactly at the Fermi energy Ep, where it is prob- 
ably a good approximation. 

This follows from considerations on the Anderson im- 
purity model. In the case of a single localized impurity 
orbital, the exact solution gives the same phase shift and 
density of states at Ep as the restricted mean-field so- 
lution (i.e., the mean- field solution where the two spin 
populations are constrained to be equal to each other, 
which is basically equivalent to a non-magnetic calcula- 
tion). In a more general case the Friedel sum rule reads 

AZ = --Y 8 lma (E F ) (2) 

IT ^-^ 

Imcr 

where AZ is the valence difference between impurity and 
host while Si m<J (Ep) is the phase shift at Ep for angular 
momentum Im and spin o and is directly related to the 
scattering matrix and scattering rate. If one assumes 
that only the d electrons screen the impurity charge 
of a transition-metal impurity and that all d orbitals 
are equivalent and have the same phase shift 5 c i(Ep), 
and additionally takes into account that Sd(Ep) is spin- 
independent in the Kondo regime, then the Friedel sum 
rule reads AZ = — (10/ it) Sd(Ep). Thus the phase shift at 
Ep is fixed to the same value whether one considers the 
restrictcd-mcan-field or the exact solution. In reality the 
assumption that only the d electrons screen the impurity 
is an approximation, as is the assumption that the phase 
shift of all d states is equivalent, since there is a small 
crystal-field spitting. In practice, however, the two as- 
sumptions are relatively well fulfilled, as one can see from 
comparison of experiment with calculations on the resid- 
ual resistivity of simple-metal dilute alloys with 3d impu- 
rities, where the impurity electronic structure was calcu- 
lated within the density-functional theory and restricted 
to be non-magnetic (see, e.g., Refs. [27l l 29l ). Therefore, the 
results obtained here by the non-magnetic calculations 
should be considered to be a reasonable approximation, 
since only Fermi-level phase shifts enter the calculation. 
The approximation clearly breaks down for properties 
that probe the phase shift off the Fermi level, e.g. in 



3 



the thermopower where the derivative dSi (E) / dE at E-p 
enters. In this case the density-functional calculations^ 
cannot account for the giant thermopower encountered 
for the Kondo impurities Mn and Fe in Al. 



II. THEORY 

Details of the KKR method for electronic structure 
calculations and impurity scattering have been published 
elsewhere.— ~— Here we focus on the definition of some 
quantities which will be useful in the discussion. 

In the KKR method, the Bloch wavefunction ipz(r + 

R n ) at a point r within an atomic cell n sited at R n is 
expanded as 



^(r + R n ;E) = Y,o %nL {E)R nL {r;E) 



(3) 



where the R n L (r; E) are regular scattering solutions of 
the Schrodinger equation of the atom n embedded in free 
space with boundary condition of an incoming spheri- 
cal wave of angular momentum L = (Z,m), while E de- 
notes the energy eigenvalue. The coefficients C£ nL (E) 
are found by the solution of the KKR secular equation 
for lattice site and propagated by a Bloch factor to any 
other site. 

In the presence of an impurity, an incoming Bloch wave 
^ evolves after scattering into an outgoing wavefunction 

■)/;™ p (with k here denoting the initial incoming state). 

The associated transition amplitude (T matrix), 



differ only by the complex conjugation of R n L(r;E)). 
Since At r l L ,{E) and A2 i ,(£J) vanish at sites where the 
perturbation of the potential, AV n (r), is zero, the sum- 
mation in Eq. © is limited to the sites for which the 
self-consistent potential is appreciably affected by the 
presence of the impurity; usually the impurity site and 
its nearest neighbors are sufficient. Only these sites are 
taken into account also for the calculation of the impurity 
structural Green function via an algebraic Dyson equa- 
tion 
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from the host Green function Q°^' nn (£;) 

Finally, the probability for an electron scattering from 
a state k to a state k' is given by 

where N is the total number of atoms in the crystal and 
c is the impurity concentration (iVc is the number of ran- 
domly positioned impurities). Assuming that each impu- 
rity scatters independently, the lifetime of a state k can 
be found by summing up the probabilities for scattering 
of a state k into all possible final states k! 

r % - E p & - t Nc T, fe I 2 - e p) ■ en) 



T & (E) := J d A r ^, (r; E)AV(r)^(r; E), (4) 
is calculated using the expression 

r tt (£)=EE c Ll ( E ) T ll' {E)c Pn , L , (E) (5) 

nn' LL' 

where the algebraic form of the T matrix 

Tit ==X>"" k^^^ + gG^r'A^^^ 

(6) 

contains information on the impurity structural Green 
function G L £,' nn (E) and on the difference in the atomic 
scattering properties between the embedded-impurity 
and substituted-host atom at site n: 



At n LL ,(E) := / d A r R nL (r; E)AV n (r)R^,(r; E) (7) 

J cell n 

Al L ,(E) := f d 3 r R* nL (r; E)AV n (f)R^,{r- E) , (8) 



with R 1 ™^/^] E) being the regular scattering solutions of 
the impurity atom in free space in analogy to the host- 
atom solutions R n L(r;E) (the r.h.s of Eqs. (JT]) and © 



The summation over k' can be transformed to an integral, 
such that the inverse lifetime is obtained by integration 
over the Fermi surface S(E-p) 



k 



1 2irN 2 c 
Vbz h 2 



f ds k' 
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ls(E F ) v k> 





(12) 



where dS^, is the Fermi surface element, is the Fermi 
velocity and Vbz the Brillouin zone volume. 



III. COMPUTATIONAL AND NUMERICAL 
ASPECTS OF THE CALCULATIONS 

The electronic structure of the investigated finite- 
thickness films and defects has been calculated self- 
consistcntly within density-functional theory using the 
principal layer technique implemented in the KKR Green 
function metho d i 35 using the experimental lattice pa- 
rameter and an angular momentum cutoff of Z max = 3. 
Exchange and correlation effects were included within the 
local-density approximation to density-functional theory 
in the parametrization of Vosko et al,— while relativistic 
effects were taken into account in the scalar relativis- 
tic approximation (ignoring spin-orbit coupling). Ex- 
cept from section UlI CI where the difference between full 
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potential (FP) calculations and those performed in the 
atomic sphere approximation (ASA) is compared, the 
ASA was used. The perturbed region, where charge re- 
laxation is allowed, is restricted in our calculations to a 
cluster of 13 sites, thus the shell of nearest neighbors, 
unless otherwise stated. Structural relaxations are not 
taken into account, as they could have a quantitative 
but not qualitative effect to the trends that are related, 
as we find, to symmetry reduction, reduced screening of 
adatoms, and the variation of localization of the surface 
states as a function of film thickness. 

Our calculations are performed on finite-thickness 
films. The situation of a half-infinite crystal with a sin- 
gle surface is approximated by increasing the film thick- 
ness up to 40 layers. Furthermore, the vacuum region 
is described by empty atomic sites at 3 (or in the FP 
calculations 4) layers above the surface; these are called 
'vacuum layers', composed of 'vacuum sites'. 



Fermi surfaces of copper, silver and gold (111) 
films 



The three investigated materials, copper, silver and 
gold crystalize in the fee structure, have a very similar 
electronic structure and therefore their Fermi surfaces are 
very similar to each other, too. As these systems are char- 
acterized by a two-dimensional periodicity, their Fermi 
surface consists of one-dimensional curves that form ring- 
like structures, except in the vicinity of the Brillouin-zone 
boundary, where hexagonal- like structures occur. The 
number of rings on the Fermi surface scales linearly with 
the number of layers, as can be seen in Fig. [TJ for the 
example of a copper (111) Fermi surface for films with 
6 (left) and 18 copper layers (right). The two innermost 
rings represent the surface bands which arc formed in the 
gap of the surface-projected band-structured 

Surface states are localized at the atomic layers close 
to the surface and decay exponentially into the bulk and 
into the vacuum. Since in a finite film there are two sur- 
faces (say "Left" and "Right"), also two surface states 
V>l and ipR appear and interact, forming a bonding and 
an anti-bonding hybrid, %j)± = (-0L ± V'r)/v / 2 with ener- 
gies E±(k) (see Fig. [2]). The coupling of these two states 
manifests in a splitting between the two inner rings of 
the Fermi surfaces, decreasing with increasing film thick- 
ness since the overlap of the two surface states decreases 
exponentially. In the sections on the lifetime calculation 
we show results on the innermost surface-state ring at 
k x > 0, k y = 0, noting that the lifetime Tg changes only 
very weakly for different k corresponding to the same 
surface-state ring, due to the isotropic character of the 
surface states. 

In the limit of infinite film thickness only one of the two 
surface bands survives at each surface (either or ^r). 
However, the overlap of these states with the surface- 
impurity is higher by a factor y/2 compared to the over- 
lap of tp± with the impurity. Therefore, for comparison 



with experiments performed not on ultrathin films but 
on crystal surfaces, the scattering rates which are calcu- 
lated here in the ip± basis can be taken from the thick 
but finite film calculations but have to be multiplied by 
a factor 2. 
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FIG. 1: Fermi surfaces of Cu(lll) films direction for 6 (left) 
and 18 layers (right) shown within the first Brillouin zone. 
In the calculations three 'vacuum layers' have been added on 
each side. The splitting of the two surface states (represented 
by the two innermost rings) decreases with the number of 
layers such that it is not visible in this scale for the 18 layers- 
film (the two innermost rings fall practically on top of each 
other and appear as a thicker circle). 




layer index 

FIG. 2: Charge (in units of the electron charge) per atom 
in the different layers for the two surface states of a copper 
film with six layers and for the innermost Fermi ring of a 
copper film with 40 layers (the three outmost layers corre- 
spond to vacuum). Here, |+) and |— ) denote the bonding 
and antibonding surface state, i.e., corresponding to the in- 
nermost/second innermost Fermi ring, respectively. For the 
thinner film, the charge in each layer is larger than for the 
40 layer-film, because the number of layers in which the sur- 
face states can penetrate is limited. The two states of the 40 
layer film are almost degenerate and nearly indistinguishable 
in their charge-per-layer picture. 
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B. Testing of the lifetimes via the optical theorem 

Contrary to the case of bulk systems where 
momentum-relaxation times and residual resistivities due 
to scattering at impurities have been calculated in ab ini- 
tio calculations (e.g. in Ref. I26ll27ll38[ ) . there are no nu- 
merical or experimental references to compare our results 
of lifetime reduction of surface states caused by adatom 
scattering. However, another possibility to attest the nu- 
merical correctness of the results is given by the optical 
theorem, according to which the elements of the scatter- 
ing matrix T^^, must obey 

h lmJ ^" v BZ h J S{E%) hv p |ifefe '' ~ T * ■ 

(13) 

This identity turns out to be a very sensitive consistency 
probe. Also, it is numerically fulfilled only if the numer- 
ical convergence of the quantities involved is very good. 
We have calculated surface-state scattering rates tC 1 for 
impurities in the first surface layer as well as in adatom 
position on top of the surface at the fee threefold-hollow 
site on a 6 layers Cu(lll) film via both sides of Eq. (fT3"|) . 
The deviation between the left and the right side of 
Eq. (fT3| was in all cases less than 1.5% (the latter devi- 
ation was obtained for surface-state scattering rates off 
impurities in adatom position in FP calculations). The 
accuracy of the calculations is very sensitive to the num- 
ber of k points used for the calculation of the structural 
Green function of the host, GYl'(E), which is used to 
calculate the impurity Green function Q™v> nn v j a a 
Dyson equation. Additionally, four vacuum layers are 
necessary to obtain the above mentioned accuracy. 



C. Comparison of ASA vs. FP and single-site 
vs. multiple-site calculations 

The computational effort in ASA is lower than in FP 
calculations. However, at surfaces, where the symme- 
try is broken, non-spherical components of the potential 
could play an important role. Therefore, for the exam- 
ple of a film of 6 layers of copper, surface-state lifetimes 
within ASA and FP calculations are calculated and com- 
pared to each other. 

The results are presented in Fig. [3] In the two pan- 
els on the left the self-consistent perturbed potential of 
the impurity and the nearest neighbors (13 sites in total, 
which we call a multiple-site calculation) has been in- 
cluded in the calculation of the scattering matrix TVv, . In 
the two panels on the right, calculations in the single-site- 
approximation are presented, i.e. only the impurity site 
is taken in the summation of Eq. ([5]) , but with the impu- 
rity potential previously converged using also the nearest- 
neighbors perturbation. All scattering rates, i.e. inverse 
lifetimes 7V~ , are given in units of ps™ 1 per atomic per- 
cent of impurities. 



For impurities in the first surface layer (the results are 
shown in the two upper panels of Fig. [3]), the agreement 
between the ASA and FP schemes is very good, both for 
single-site and multiple-site calculations. This is not the 
case for adatoms as can be observed in the two lower pan- 
els of the same figure where a quantitative difference of 
maximal 25% can be observed. Nevertheless, the trends 
in FP and ASA are qualitatively similar. Since we are 
interested in these trends, we will restrict our further cal- 
culations to the ASA. 

The qualitative difference between single-site multiple- 
site calculations in the case of adatoms raises the ques- 
tion whether the cluster consisting of nearest neighbors 
is large enough to obtain converged results. Therefore, 
multiple-site calculations (not shown here) including ad- 
ditionally the shells up to the third nearest neighboring 
sites (43 atoms in total) have been performed. They show 
that the results remain practically stable beyond the 
nearest-neighbor approximation, which we adopt hence- 
forth. 



IV. ANALYSIS OF SURFACE-STATE 
LIFETIMES OF Cu FILMS 

The most striking feature of our results is the qualita- 
tive difference between the trends of scattering rates at 
adatoms and those of impurities in the surface, as can 
be seen in Fig. 03 While the trend for impurities in the 
first surface layer (lower panel) resembles that which has 
been found for scattering at impurities in bulk , 26 ' 27 the 
trend for the scattering rate at adatoms (upper panel) 
does not. 

In the case of impurities in the first surface layer as well 
as in bulk, the observed trend follows the local density of 
states of the impurities (see Fig.@|: For the d scatterers, 
the scattering rate shows a peak with a maximum at 
V and Cr as the d resonance crosses the Fermi energy. 
Then, as the p states start crossing Ep, the scattering 
rate increases again starting off quadratically according 
to Linde's rule . 39 ' 40 

Basically, this correlation between the impurity local 
density of states and the scattering strength of impurities 
can be traced back to conditions which are not necessarily 
met for adatoms. The first condition is that the screen- 
ing of the impurity charge in metals takes place mostly 
in the impurity atomic cell. Thus the perturbation of the 
potential in nearby cells is rather small, while the impu- 
rity wavefunctions are relatively atomic-like in the cell 
and matched to phase-shifted host waves outside. The 
second condition is that the free-electron-likc behavior 
in the bulk of normal metals together with the nearly 
spherical symmetry of the impurity potential allows for 
independent scattering of the s, p and d without mutual 
interference. Thus the angular-momentum resolved lo- 
cal density of states reflects the scattering properties of 
the angular-momentum channels independently to each 
other. 
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FIG. 3: Scattering rates tC 1 for impurities in the first surface layer (upper panels) and adatoms (lower panels) obtained within 
the atomic sphere approximation and full-potential calculations. While the agreement between the two calculations is very good 
for impurities in the first surface layer, for adatoms the two calculation schemes lead to qualitatively similar but quantitatively 
different results. In the left panels results of multiple-site calculations are shown, i.e., with a cluster of the nearest neighbors 
included in the potential perturbation. In the right panelsthe comparison is shown for single-site calculations. 



For adatoms, however, the screening is not as efficient, 
thus the atomic size can increase significantly. Here, by 
atomic size we mean the volume in which the impurity 
induces a significant potential perturbation. In addition, 
the reduced symmetry allows for interference between dif- 
ferent angular momenta, notably (as we will sec) between 
s and p z - Impurities in the surface layer are apparently 
close enough to the impurity-in-bulk situation that they 
behave in the same way. 



Having in mind these comments, we proceed to an 
analysis of the results. 



A. Adatom scattering 

Considering the trend of adatom-scattering rates in 
the multiple-site approach (see the left-bottom panel of 
Fig. E}, there are two features to observe: The first part 
of the trend (from Sc to Fe/Co) is dominated by a rel- 
atively large scattering rate, showing a plateau which is 
absent in the single-site calculations (right-bottom panel 
of Fig. [3]). On the other hand, the second part of the 
curve (from Ni to Se) is qualitatively the same for single- 
site and multiple-site calculations. However, it behaves 
in an uncommon manner, showing a kink-like structure, 
which is completely different from the trend for the sp 
impurities in the surface. 

The plateau in the scattering rates of the early 3d 
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FIG. 4: Top: Local density of states of Sc, Cr, Ni and Ga 
adatoms on a Cu(lll) film of 40 layers. The density of states 
does not differ qualitatively from that of impurities in the 
surface layer. Bottom: Density of states summed over the 
vacuum sites surrounding the adatoms. 



scatterers can be understood by the significant differ- 
ence between the single-site approximation compared to 
multiple-site calculations (see Fig. [SJ, which hints at an 
atomic size effect. These atoms (Sc, Ti, V) have rel- 
atively large atomic radii (which are r at = 1.62 A for 
Sc, 1.45 A for Ti, 1.34 A for V)ii and therefore extend 
much more into the vacuum than e.g. a Ni adatom with 
r at = 1.24 A. This is not any more the case for atoms em- 
bedded in the surface or in bulk, as the strong screening 
strongly reduces the atom size. In order to verify this line 
of arguments we have performed additional calculations 
where - apart from the adatom site - only the perturbed 
surrounding vacuum sites or only the perturbed substrate 
Cu atoms are included in the expression for T^p (see the 
summation over sites n, n' in Eq. (O). These calculations 
reveal that the largest contribution next to the adatom- 
site arises from the surrounding vacuum potential, and 
not from the Cu atoms in the surface layer. We also ex- 



70 
60 
50 
40 
30 
20 
10 

-10 
-20 



incl. nearest neighbors — \— 
single site — --X-- 
rence 



diffl 



.-X 



X 



.X 



Sc Ti V Cr Mn Fe Co Ni Cu Zn Ga Ge As Se 
adatom 



Sc 






Cr 






Ni 






Ga 






unperturbed 






















< 










Sc Ti V Cr Mn Fe Co Ni Cu Zn Ga Ge As Se 
adatom 

FIG. 5: Top: Surface-state scattering rate r^ 1 in ps _1 
per atomic percent for scattering at adatoms. Red curve: 
Multiple-site calculations. Green curve: Single-site approxi- 
mation. Blue curve: Difference between the two cases. Bot- 
tom: Excess charge (in electron-charge units) induced by 
adatoms in the surrounding vacuum sites. For the d adatoms, 
the trend observed in the difference between the scattering 
rates calculated including nearest neighbors and single-site 
scattering is similar to the trend of the charge induced by the 
adatoms. 



amined the density of states at E-p around the impurity, 
shown in the lower panel of Fig. [4] We found that, in- 
deed, there is an increased DOS in the vacuum around 
the early 3d adatoms, while this contribution is reduced 
for the late 3d series. When the p states start crossing E-p 
(at Ga) the size of the atom grows again, but apparently 
the qualitative trend is governed by different effects (sec 
next paragraph), although there are some quantitative 
differences between single- and multiple-site calculations. 
The larger extent of the early 3c? adatoms also manifests 
in a larger charge transfer to the surrounding vacuum 
sites. This becomes visible in the lower panel of Fig. [SJ 
where we show the charge AQ in the surrounding vacuum 
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FIG. 6: Angular-momentum resolved scattering rates for 
adatoms on a 40 layer Cu(lll) film. The inclusion of just 
the s and the p z channel in the single-site scattering matrix 
Tll' (compare Eq. ^) for the 4sp adatoms is already a good 
approximation for the total scattering rate. Scattering in only 
the s channel by far overestimates the total scattering rate. 
The same can be observed for Ge, As and Se adatoms when 
scattering is restricted to the p channel. The total scatter- 
ing rate is then reduced due to a destructive interference of 
s and p scattering, which can be verified by the investigation 
of Friedel oscillations (see Fig. . 
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sites induced by the adatom. For the 3d adatoms we find 
a decreasing trend for AQ similar to the difference be- 
tween the scattering rates including the nearest neighbors 
and single-site calculations (compare the dotted line of 
the upper panel with the lower panel of Fig. [5]). We note 
that, in the spirit of the self-consistent density-functional 
approach, the excess charge causes extra scattering only 
through the induced change in the potential. 

In order to understand the second part of the trend for 
scattering at adatoms (from Ni up to Se), which shows a 
zig-zag structure, we analyze the single-site results, be- 
cause qualitatively they do not differ from multiple-site 
results, while they allow for an angular-momentum re- 
solved investigation of scattering rates. In the spirit of 
such an /m-decomposition, we restrict the single-site con- 
tribution n = n' = of the scattering matrix T^£, de- 
fined in Eq. ([6]) to different L channels. The results are 
presented in Fig. [5] for T s ° s ° and together with the 
total single-site contribution. Due to the presence of the 
surface, which breaks the symmetry, we expect that inter- 
ference is possible among s and p z waves (the z direction 
is taken along the surface normal) . The inclusion of the ,s 
and the p z channel already leads to a curve very close to 
the total single-site result (d z i waves are insignificant for 
sp adatoms). However, restriction of T 00 to the s channel 
only overestimates the scattering rate for adatoms from 
Zn to As, while restriction to the p channel only overes- 
timates the scattering rate for Ge, As and Se adatoms. 



FIG. 7: Friedel oscillations in the vacuum near the surface and 
as a function of the distance from the impurity calculated for 
Ga and Ge adatoms. The difference of the local density of 
states as a function of distance to the adatom site is calcu- 
lated, i.e. n^P^F)-^ 034 ^), with nl ost (E F ) the density of 
states of the vacuum site of the host and n™ p (i?F) the density 
of states of the same site perturbed by the adatom. Here, the 
index L means that the scattering has been artificially con- 
fined to the L channel only (s or p z ). Both for Ga and Ge 
adatoms a destructive interference of the s and p z channel is 
observed as they are out of phase, d scattering is found to be 
very small and can be neglected. 



The conclusion is that there is destructive interference 
between the s and the p z channels which is allowed by 
the breaking of symmetry due to the surface. The zig- 
zag structure, with a peak at Zn and a dip at Ge, occurs 
because the maxima of the s and p channels arc phase- 
shifted with respect to each other. 

For a more elaborate analysis of interference effects we 
have calculated Friedel oscillations. In particular we look 
at the difference of the local density of states An^Ep) 
in the vaccum close the impurity, and in a direction par- 
allel to the surface, between the system with impurity 
and the host system. To this end we integrated, for the 
two cases, the density at Ep over the atomic spheres of 
the vacuum sites with and without the impurity, obtain- 
ing An L (£: F ) = n h L ost (E F ) - n^ p (E F ), for Ga and Ge 
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adatoms. The index L here means that we made again 
an angular-momentum decomposition, by allowing scat- 
tering by only a certain L channel while setting the oth- 
ers to zero. In this way we are able to see, in real space, 
the relative phase shift in channels which are expected 
to interfere, e.g., between the s and p z channel. This 
analysis is presented in Fig. [7] as function of the distance 
from the adatom site. We see that oscillations in the s 
and p z channel are mutually phase-shifted and interfere 
destructively. 



B. Separation of scattering into bulk and surface 
states 

After a scattering event, an electron that was previ- 
ously occupying a surface state fc S urf can end up in an- 
other surface state k' sur f or in a bulk state £; b uik- The 
distinction between the two cases can be of importance 
to effects such as surface-state mediated interactions be- 
tween defects, surface resistivity, the lifetime of surface 
plasmons, etc.. The method used to calculate the surface- 
state lifetimes given in section IPTl Eq. (|T2"j) . allows to dis- 
tinguish between the two cases. The total scattering rate 
of a (surface) state characterized by a wavevector k is 
composed of a contribution 1 / t~ given by 



r surf 
k 



ds % , 



S(E F ) 



± kk' I 



(14) 



where the integration is performed only over the surface 
bands and analogously a contribution l/rS ulk , where only 
the bulk end-states of the Fermi surface are taken into 
account: 



1 



r bulk 

k 



2ttN 2 c 

v BZ h 2 



S(E F 



r kk' 



(15) 



Obviously, considering the definitions of tv, t5 U1 ' and 



5 the relation 



t~_ -surf 



r bulk 



(16) 



holds. 

The two contributions have been calculated for films 
consisting of 6 and 40 layers of copper for impurities in 
the surface layer as well as for adatoms and are presented 
in Fig. [H The most important observation is that 1 /T^ ulk 
and 1 / T"5 are on the same order of magnitude, even in 
the case of adatoms. Naturally, variations occur as a 
function of film thickness, which we will now address. 

Qualitatively, we expect that the relative scattering 
rate into bulk and surface states should depend on (i) the 
amplitude of the surface states, ips(Rimp), at the impu- 
rity site, (ii) the amplitude of the bulk states, "0B(i?i mp ), 
at the impurity site, and (iii) the available number of 



bulk states, Ab- About point (i), as a function of the 
number of layers in the film, Al, ips (-Rimp) first de- 
creases before reaching saturation, since in very thin films 
the number of layers in which the surface states can 
penetrate is limited (see also Fig. [2]). From this effect 
alone we expect an overall reduction of both 1 /t~ ui1 ~ 

|(^(i?imp)|T|^(i?imp))| 2 



IV-s(#im P )| 4 and l/rS ulk 



-imp;^ irin-'-imp;/! iV's(-Rimp) 1 2 with increas- 
ing TVl before a saturation takes place, but with l/r| urf 

being more strongly affected due to the 4-th power of the 
amplitude. About point (ii), "0s(i?i mp ) ~ 1/\/Al, yield- 
ing a factor of 1/Al to l/rS" 1115 ; however, this reduction 

is partly compensated by point (iii), which contributes a 
factor A B = A L - 1 to l/rS uIk . 

' k 

Combining (i), (ii) and (iii) gives a behavior of the 
form l/rS ulk ~ |i/; s (i? imp )| 2 (l -1/N L ) for scattering into 
bulk states, which should show a saturation for large Al, 
but could be either increasing or decreasing at small Al 
depending on the relative importance of the two terms. 
On the other hand, point (i) gives a behavior of the form 
l/r| ~ IV'S^imp)! 4 for scattering into surface states, 
which should be first decreasing with Al and then also 
saturating. 

The qualitatively expected behavior is consistent with 
our calculations (Fig. [8]). The scattering rates into sur- 
face states are larger for the 6 layer film than for the 40 
layer film by approximately a factor 2, while the scat- 
tering rates into bulk states are a little larger for the 40 
layer film. 



COMPARISON BETWEEN Cu, Ag AND Au 
FILMS 



In the previous section we have analyzed surface-state 
lifetimes for impurities in and on copper films. Here we 
extend our investigations to silver and gold films. Even 
though the band structure of the three noble metals is 
similar, differences are anticipated mainly due to the 
change of energetic position of the surface state, espe- 
cially in Ag. 

The calculated data, shown in Fig. [9j are qualitatively 
similar for all three hosts. Scattering rates off adatoms 
are largest for the first elements of the row and remain 
almost constant until Mn; just as in the case of the cop- 
per film (see previous section), this should be due to the 
larger atomic radius entailing a larger extent of these 
adatoms into the vacuum. Furthermore, all three host 
materials show a clear trend that scattering at adatoms 
on a film of six layers is enhanced compared to that of 
larger thickness. The reason is the higher localization 
of the surface state for thin films, as explained in the 
Cu(lll) case. 

We show here results for the innermost surface state 
only, whose lifetime, especially for the 6 layer films and 
for some adatoms, remarkably differs from that of the 
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FIG. 8: Scattering rate 1 for the innermost surface state together with the contributions of scattering rates to bulk and 

surface states l/rS ulk and l/r2 urf for films with 6 layers (left panels) and 40 layers of copper (right panels) for adatoms (top) 
and impurities in the first surface layer (bottom). 



second surface state. While for the early 3d elements 
the difference is relatively small (a few percent), it can 
increase up to 60% for Ge and 48% for As adatoms. For 
the 40 layer film, this difference almost vanishes as does 
the splitting between the two surface bands. 

Concerning scattering at impurities in the surface, the 
situation is somewhat complicated. Qualitatively, for all 
three host materials a clear maximal scattering rate for 
the 3d elements with half-filled shells (V, Cr) is observed 
as well as high scattering rates for the sp scatterers. 
Hence, as expected, the global trend reflects the situation 
of scattering at impurities in the bulk. However, large 
quantitative differences among the three host materials 
are observed when considering the thickness-dependence 
as well as the comparative scattering strength of 3d and 
Asp impurities. While for a copper host, the scattering 
rate does not depend much on the film thickness, espe- 



cially for silver films the thickness can change the scat- 
tering rate up to almost an order of magnitude. Actu- 
ally, such a strong variation in the Ag films is reasonable, 
because of the relatively shallow position of the surface 
state with respect to the Fermi energy.— As a result, in 
thin films quantum-confinement effects push one of the 
two surface states (the innermost ring) above Ep and the 
other lower, resulting in a significant change of the avail- 
able phase space for scattering. In thicker films the two 
states are approximately degenerate and below Ep. In 
any case, in Ag the scattering rate into surface states is 
small due to the smallness of the corresponding Fermi 
rings. 

A comparison to other theoretical or experimental re- 
sults is not possible because of lack of data; although 
these surface states have been subject of many experi- 
ments, to our knowledge no experiments have been per- 
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FIG. 9: Surface-state scattering rate Tg 1 for adatoms (left panels) and impurities in the first surface layer (right panels) for 
Cu, Ag and Au films with different number of layers. For silver films the scattering rate is most sensitive to the film thickness 
due to the fact that one of the surface bands lies above Ef in the 6 and 10 layer cases. Qualitatively the trend as a function 
of the impurity atomic number is similar for the three hosts. 



formed in which the surface-state lifetimes due to scat- 
tering at the specific impurities at the Fermi level have 
been measured. However, the order of magnitude of 
the calculated scattering rates should allow for an ex- 
perimental detection, which e.g. in ARPES reaches a 



few mcV. To compare, in inverse photoemission spec- 
troscopy a linewidth of about 23 meV ~ 35 ps _1 ?i for 
Cu, 6 mcV w 9 ps -1 fo for Ag and 21 meV ~ 33 ps~ 1 h 
for Au for electron-electron scattering processes has been 
measured^ A measurement for 1% of defects should be, 
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therefore, within the experimental resolution. 

VI. SCATTERING AT MAGNETIC 
IMPURITIES 

So far, only scattering at non-magnetic impurities and 
adatoms has been considered, in an approximation that 
was explained in the introduction. However, above the 
Kondo temperature, some of the 3d impurities become 
magnetic and scattering in the two spin channels has 
to be treated separately. We investigate now the con- 
sequences of spin polarization, showing results for the 
Cu(lll) surface only, which suffices to demonstrate the 
main effect of a double-peak structure in the scattering 
rate. 

The calculated lifetimes are shown in Fig. [10] for a 
40-laycr film. In the adatom position the 3d elements 
starting between Ti and Ni are magnetic, while as im- 
purities in the surface Ti and Ni are paramagnetic. The 
Asp elements are in all cases non-magnetic. The mag- 
netism of the 3c? elements leads to an expected double- 
peak structure in the trend of the scattering rates, which 
is already known e.g. for residual resistivity in bulk above 
the Kondo temperature^— and originates from the off- 
set of the d resonances of the two spin channels, which are 
mutually repelled by the exchange interaction. Thus, a 
first peak of the scattering rate is observed for Ti, where 
the d resonance of the majority channel is centered at the 
Fermi level, while a second peak appears for Fe and Co 
in the minority-spin channel. For Cr impurities, where 
scattering rates are large in the case of paramagnetic im- 
purities, the scattering rate is most drastically reduced 
as E-p lies between the resonances. 

VII. SUMMARY 

In summary, we have presented a systematic study of 
the lifetime of the (111) surface state of the noble metals 



considering scattering at the Fermi energy due to 3d and 
Asp impurities atoms on and in the surface. 

Our main finding is a qualitatively different trend, as 
a function of the impurity atomic number, for adatoms 
compared to surface-embedded impurities. While the lat- 
ter case is similar to the trend for impurities in the bulk, 
reflecting the sharpness of the d-resonance and Lindc's 
rule as the d and p states cross the Fermi level, adatom 
scattering is found to be dominated by other criteria such 
as the size of the adatom and interference effects between 
the scattering amplitudes of different angular momenta. 
The conclusion is that the surface strongly affects the 
scattering properties by the reduction in charge screen- 
ing as well as by the lowering of symmetry. 

A separation of scattering rates into surface-to-surface 
and surface-to-bulk parts shows that the two contribu- 
tions are comparable in magnitude for adatoms as well 
as for embedded impurities. Additionally, in the case 
of ultrathin films we have found that the scattering into 
surface-states contribution decreases with film thickness, 
since the surface state extends further into the bulk in 
thicker films, and finally saturates. Furthermore, we have 
found that the lifetime is sensitive on the film thickness, 
especially for Ag. This is due to the very shallow sur- 
face state of Ag, which can be shifted in energy due 
to quantum confinement effects, changing the available 
phase space for scattering. 

The calculated lifetimes should be detectable as a 
surface-state linewidth within the resolution of photoc- 
mission spectroscopy already for impurity concentrations 
1%. 
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